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Abstract 

We study the presence of a logarithmic time scale in discrete approximations of 
Sawtooth Maps on the 2-torus. The techniques used are suggested by quantum 
mechanical similarities, and are based on a particular class of states on the torus, 
that fulfill dynamical localization properties typical of quantum Coherent States. 
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1. Introduction 

Under the term of Quantum Chaos goes a rich phenomenology of behaviours [1-3] proper 
to quantum systems whose classical limit presents typical chaotic features as positive 
Lyapunov exponents (hyperbolic regime) [4-6]. 

The footprints of classical chaos are usually studied semi-classically when a suitable 
" like quantization parameter goes to zero; one then examines the differences between 
quantum and classical behaviours. In the hyperbolic case, quantum chaos reveals itself 
through the presence of a time-scale, over which quantum and classical motions mimic 
each other, that increases as — log"/i" [1-3,7-9]. This peculiar logarithmic time scale has 
to be compared with the scaling "h"~ a , a > 0, which is proper of quantum systems with 
regular classical limit [1]. 

Heuristical explanations of the logarithmic time-scale already indicate that the phe- 
nomenon is not exclusive of quantum systems, and thus of non-commutativity, but that 
it should also be present when the classical dynamics is looked at as the continuous limit 
of a family of discrete classical systems. [10]. 

Intrinsically discrete systems [11] and discretized classical continuous systems [12- 
14] have recently been objects of numerical analysis concerning the entropy production 
and the presence of a logarithmic time scale, whereas the ergodic properties of discretized 
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discontinuous maps have been addressed in [15]. 

In the following, we shall rigorously show this fact to be true for Sawtooth Maps on 
the 2-dimensional torus [16-18]: this will be done by forcing them to move on a square 
lattice and by retrieving the continuous dynamics when the lattice spacing goes to zero. 
Because of the analogies between quantization and discretization, we will make use of 
technologies strictly resembling the so-called Anti-Wick quantization [19]. 

We shall prove that a time-scale logarithmic in the lattice-spacing appears; in com- 
parison to previous results obtained studying numerically the entropy production [14], 
a rigorous continuous limit is established that succeeds in controlling the discontinuities 
of Sawtooth Maps. Despite their classical nature, the entropy previously investigated 
was quantum mechanical; somewhat analogously, in this article, Sawtooth Maps will be 
studied by means of states, which play a role similar to quantum Coherent States, whose 
choice is naturally provided by the lattice structure of discretized Sawtooth Maps. They 
will be shown to satisfy a dynamical localization property that makes them remain local- 
ized around the trajectories of the continuous dynamics, but only on a logarithmic time 
scale. 

2. Classical Dynamical Systems 

Classical dynamics is usually described by means of a measure space X ', the phase-space, 
endowed with the Borel cr-algebra and a normalized measure /i, n{X) = 1. The "volumes" 



of measurable subsets E C X represent the probabilities that a phase-point x S X belong 
to them: the measure [i defines the statistical properties of the system and represents 
a possible state, which is taken to be an equilibrium state with respect to the given 
dynamics. 

In such a scheme, a reversible discrete time dynamics amounts to an invertible 
measurable map S : X ^ X such that fi o S = fi and to its iterates {S k \ k G Z}: 
phase-trajectories passing through x S X at time are then sequences {S k x} [6]. 

Classical dynamical systems are thus conveniently described by triplets (X,/j,,S); 
in the present work, we shall focus upon the following choices: 




X: the 2-dimensional torus T 2 = R 2 /1? = {x = {x 1 ,x 2 ) G K 2 (mod 1)}; 
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//: the Lebesgue measure, fj,(dx) = dxidx2, on T 2 ; 

5: an invertible measurable transformations on T 2 that preserves the Lebesgue mea- 
sure. 

It is convenient to associate an algebraic triple to the measure-theoretic 

triple (T 2 ,/i, S), consisting of 

M: the (Abelian) Von Neumann *-algebra (T 2 ) of essentially bounded functions 
on T 2 [20,21]. 

oj^: the state (expectation) on M, given by 

Wm : L™ (T 2 ) 3 f ,— > u^f) := / fi(dx) f(x) G R+ • (1) 

Jt 2 

0: the automorphism of M. such that (/) = /oS,wo6=w. 

In the following, we shall consider a discretized version of (T 2 ,//, S) which arises by 
forcing the continuous classical system to live on a square lattice Ln C T 2 of spacing jj: 

Ln '- = {n \ p e ( Z / NZ ) 2 } > ( 2 ) 
where (Z/iVZ) denotes the residual class (mod A), that is ^ pi ^ N — 1. 

Taking the A 2 points as labels of the elements {| i^^/nt) 2 of an orthonormal 
basis (o.n.b.) of the N dimensional Hilbert space Htf, ftf:= A 2 , we will consider discrete 
algebraic triples \T>Af,Tj\f,@j\f), consisting of 

Vj\f. an M x M matrix algebra diagonal in the orthonormal basis introduced above; 
Tj\f. the uniform state (expectation) on T>n defined by 

t M :V M B D^T M (D):=^Tr(D)GR + ; (3) 

Qj\f. an automorphism of Vj^f suitably reproducing when A — ► oo (see Section 4.2). 
Remark 2.1 

As it will become evident in the following, up to a certain extent, discretization 
resembles quantization; in the latter case, instead of T>_\f, one deals with non- 
commutative matrix algebras, the typical instance being the finite dimensional 
quantization of the Arnold Cat Map [22,23]. 
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3. Discretization of phase-space 

As sketched in the previous Remark, we proceed now to setup a discretization procedure 
close to the so-called Anti-Wick quantization [19]. 

Given the classical algebraic triple {U^ (T 2 ) ,0^,0), the aim of a discretization- 
dediscretization procedure (specifically an AA-dimensional discretization) is twofold: 

• finding a pair of *-morphisms, Jx,oo mapping L£j° (T 2 ) into the abelian finite di- 
mensional algebra T>x and Joo,N mapping backward £y into L™ (T 2 ); 

• providing an automorphism Qjy, the discrete dynamics, acting on Vj\f such that it 
approximates the continuous one, 0, on (T 2 ) as follows 

J^oQ^oJ^^Qi • (4) 

The latter requirement can be seen as a modification of the so called Egorov's property 
(see [24]). Intuitively, a discrete description of the measure-theoretic triple (T 2 ,/x, S) 
becomes finer when we increase N, the number of points per linear dimension on the grid 
Ltv in (2): this corresponds to enlarging the dimension of the Hilbert space Tlj\f associate 
to the corresponding algebraic triple (T>n, t/v, Qj\f) . In this sense, the lattice spacing 
a:= jj of the grid Ln is a natural "discretization parameter" playing an analogous role 
to the quantization parameter h. 

The difficulty is to find convenient *-morphisms JV.oo and Joo,Af that set up a 
rigorous asymptotic (in N) correspondence, of functions on (T 2 ) and matrices in T>n 
and, above all, between the discrete dynamics Qj\f and the continuous one 0. 

Due to the similarities with quantization, we shall consider a discretization procedure 
based on states that we shall refer to as Lattice States (LS for short) which mimic the use 
of Coherent States in the study of the semi-classical limit. In the next section we will 
give a suitable definitions of LS belonging to the Hilbert space TCj\f, that we shall use to 
discretize (L~ (T 2 ) ,w m ,0). 

3.1. Lattice States on T 2 

In analogy with the the properties of quantum Coherent States, we shall look for a class 
{\Cj\f(x)) | x £ T 2 } £ TCj\f of vectors, indexed by points x £ T 2 , satisfying the following 
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conditions which are borrowed from analogous quantum ones [25]: 
Properties 3.1 

1. Measurability: x \— > \Cj\f(x)} is measurable on T 2 ; 

2. Normalization: ||CV(a:)|| 2 = l,x€ T 2 ; 

3. Completeness: J\f / fj,(dx) \Cjt/(x))(Cj^(x)\ = 1; 

Jt 2 

4. Localization: given e > and do > 0, there exists iVn(e, do) such that for 
N > No(e,do) and d T 2(x,y) > do one has 

M\(C M (x),C M (y))\ 2 <e. 

The symbol d T 2(x, y) used in the localization property stands for the length of the shorter 
segment connecting the two points x,y £ T 2 , namely 

Definition 3.1 

We shall denote by d V 2 (x, y) := min \\x — y + n|| R2 the distance on T 2 . 

We shall now construct a family of \Cj\f(x)}. Let [ - J denote the integer part of a 
real number, namely x — 1 < [^J ^ % is the largest integer smaller than x; further, let 
(•) denote the fractional parts, that is (x) := x — [x\ . Thus we will write 

2 ( \Nxx\ YNx 2 \ \ f (N Xl ) (Nx 2 ) \ 

3 { N ' N ) + { N ' N ) ' 

\Nx\ (Nx) 9 
or, more compactly, x = — — 1 — — . We proceed by associating to points of T 

specific lattice points. 

Definition 3.2 (Lattice States) 

Given x £ T 2 , we shall denote by xjq the element of (Z/A^Z) 2 given by 

x N = (x Njl ,x Ni2 ) ■= ( [Nxi + |J , [Nx 2 + |J ) , (5) 
and call Lattice States on T 2 the vectors \Cj\f(x)) defined by 

T 2 3 x^ \C N (x)):= \x N ) G Hm • (6) 
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Remark 3.1 



The family of states \Cj^(x)} is constructed by choosing, for each x G T 2 , 
that element of the basis of Htf which is labeled by the closest element of Ljy 
to x. 
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Figure 1: The above picture represents a square lattice (£5) of spacing | by circles and 
connecting lines. All points in the blue square I^a l) : = [ig> 10O x [to'k) c ^ 2 are 
associated with the grid point (|, |) (black dot). Thus, for all x G I/g 3\, it turns out 
that |CV(a;)) = | (3,3)) G Hat- 



Proposition 3.1 



The family of LS {|CV(a;)} satisfies Properties 3.1. 



Proof: 

Measurability and normalization are straightforward. 
Completeness can be expressed as 

N I fi(dx) {£ \C M {x)){C M {x)\ m) = tfW W,m G {Z/NZ) 2 , 

where we have introduced the periodic Kronecker delta, that is 6^q = 1 if and only if 
n = (mod iV). This is proved as follows: 

M [ fi(dx) (£\Cm(x))(Cm(x) \m) = J\f [ dxi [ dx 2 {£ \ x N )(x N \ m) = 
Jt 2 Jo Jo 
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Localization comes as follows: from Definition 3.2 (see Remark 3.1 and Figure 1), 
it turns out that \Cj\f(x)) is orthogonal to every basis element labeled by a point of 
Ltv whose toral distance d T 2 (see Definition (3.1)) from x is greater than j^j^- As a 
consequence, the quantity (Cj\f(x),Cj\f(y)) = if the distance on the torus between x 
and y is greater than Thus, given do > 0, it is sufficient to choose iVo(e, do) > \/2/do, 
to have 

N>N (e,d ) =^M {C N {x),C N {y)) = • ■ 



Remarks 3.2 



(1) The last result in the previous Proposition amounts to an even stronger 
localization property than property 3.1.4; this is due to our particular choice 
of Lattice States, which, as we shall see, is suited to the task of controlling 
Sawtooth Maps. In general, one can hardly hope to achieve orthogonality 
and must be content with the weaker localization condition 3.1.4. 

(2) Although the set of LS of Definition 3.2 fulfill Properties 3.1, which are 
typical of Coherent States, LS differ from them in that the context we are 
considering is commutative. In spite of this, it is convenient to adopt the 
formalism of Quantum Mechanics; in particular the set of LS is interpreted 
as a Hilbert orthonormal basis of Dirac kets, whose corresponding projectors 
form a partition of unit into indicator functions having support on small 
squares of the torus, as in Figure 1, whose sides scales as • 



3.2. Anti— Wick Discretization and its continuous limit on T 2 



In order to study the continuous limit and, more generally, the quasi-continuous be- 
haviour of (Vj\f, rj\f, @j\f) when iV — > oo, we follow the semi-classical technique known as 
Anti-Wick quantization. The other standard quantization technique, namely the Weyl 
procedure, despite being more straightforward and less technically heavy, is nevertheless 
more suited to smooth spaces of functions and was indeed instrumental in the study of 
discretized Cat Maps [14]. Instead, in our case, the Anti-Wick procedure is a better 
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choice due to the discontinuous character of the dynamics, as it will clearly appear in the 
next Section. 

We start choosing concrete discretization/de-discretization *-morphisms. 



Definitions 3.3 



Given the family {\Cj\f(x))} of Lattice States in 7ij\f, the Anti- Wick-like dis- 
cretization scheme (AW, for short) will be described by a one parameter family 
of (completely) positive unital map J^j - )00 : (T 2 ) — > Vj\f 

L™(T 2 )Bf^N[ f,(dx)f(x)\C M (x))(C^(x)\=:J M>00 (f)eV M . 
Jt 2 

The corresponding de-discretization operation will be described by the (com- 
pletely) positive unital map Joa^T ■ T^N ~^ Lp° (IT 2 ) 

V M 3 X ^ (C N {x),XC N {x)) =: Joosr(X)(x) G (T 2 ) . 



Remarks 3.3 



Both maps are identity preserving (unital) because of the conditions satis- 
fied by the family of Lattice States and are completely positive, since both 
(T 2 ) and P/v" are commutative algebras. One can also check that: 

\\JooM MM\\ 00 <\\g\\ 00 , geL™(T 2 ) ■ 



ii. Definition 3.3 yields t/v" ° Jm,oo — w^, with Tj\f given in (3). 



In Appendix A, more operative details are presented, whereas in the following we prove 
some simple properties that incorporate minimal requests for rigorously defining the sense 
in which the discrete dynamical systems (T>j\f, r/y, @j\f) tends to (T 2 ) , Up, O), when 



Proposition 3.2 



(1) For all / (T 2 ) and X £ V N , 

Up (g Joo,N (X)) = Tx{Jx y00 (g)*X) ; 
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(2) Forall/,<7GL-(T 2 ) 



Km TM(JM >00 (f)*JM t00 (g)) =uj /J ,(fg) = / fi(dx) f{x)g(x). 



(3) For all X G £y, and for all N G 1N + , 

JaT,oo ° Joo,N (X) = X ; 



(4) For all / G (T 2 ) 



lim j7oo,Ar ° JN,°c{f) = f M _ a.e. 
iV— >oo 



Proof: 

The first two statements in the above Proposition directly follow from Definitions 3.3 
together with (6); the latter two are equivalent and their proof can be found in [25], the 
only difference being the dimension M of the Hilbert space Hn, here M = N 2 , there 

M = n. m 



Remark 3.4 



Properties 1 and 2 in the previous Proposition show how (GNS) scalar prod- 
ucts in the discrete, respectively continuous limit, are related; properties 3 
and 4 concern instead the direct-inverse relations between the discretization 
and the de-discretization maps. 
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4. Discretization of the Dynamics 



4.1. Classical description of Sawtooth Maps 

We shall now focus on a special class of automorphisms of the torus, namely the Sawtooth 
Maps [16,17] (SM for short), that is on triples (T 2 ,//, S a ) where 




(mod 1) , a G E (7) 



1 + a i\ ( Xl y 

a l) \X2 J 

_ ( ((1 + a) ( Xl ) + x 2 )\ 
\ (a(xi)+x 2 ) J 

Remarks 4.1 

i. In the following, a point x of the torus, will correspond to an equivalence 
class of K 2 points whose coordinates differ by integer values; 

ii. without the fractional part, (7) is not well defined on T 2 for not-integer a; 
indeed, the same point a; = :E + nGT 2 ,nGZ 2 , would have (in general) 
S a (x) / S a (x + n). Of course, (•) is not necessary when a£Z; 

iii. the Lebesgue measure on T 2 is invariant for all aGl; 

iv. if a Z, the S a are known as Sawtooth Maps; 

v. when a G Z, we shall write T a instead of S a . T\ = ( \ \ ) is the Arnold Cat 
Map [6]. In general, T x G {T a } aeZ C SL 2 (Z) C GL 2 (Z) C M 2 (Z) where 
M 2 (Z) is the subset of 2 x 2 matrices with integer entries, GL 2 (Z) the 
subset of invertible matrices and SL 2 (Z) the subset of matrices with deter- 
minant one: the dynamics generated by T a G SL 2 (Z) is called Unimodular 
Group [6] (UMG for short); 

vi. after identifying x with canonical coordinates (q,p) and imposing the 
(mod 1) condition on both of them, the above dynamics reads 

q' = q + p' 

(mod 1) • (8) 

p' = p + a(q) 

This is nothing but the Chirikov Standard Map [3] in which — ^ sin(2-7r(/) 
is replaced by (q). The dynamics in (8) can also be thought of as generated 
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by the (singular) Hamiltonian 

where 6 p (t) is the periodic Dirac delta which makes the potential act 
through periodic kicks with period 1 [26]; 

vii. Sawtooth Maps are invertible and the inverse is given by the expression 

")( to -:)(:)> (9 » 

I (xi- x 2 ) \ 
\((x 2 ) - a(x\ -x 2 ))J 

or, in other words, 

q = q '~ P ' I A1\ 
(mod 1) . 

p = —a q + p' 

It can indeed be checked that S a (x)) = S' 1 (S a (x)) = x, \/x G T 2 . 
Further, 5" 1 preserves the Lebesgue measure on T 2 . 

We now list a set of properties [16-18] of Sawtooth Maps that will be used in the following 
Properties 4.1 (of Sawtooth Maps) 

(1) Sawtooth Maps {S a } are discontinuous on the subset 

7o: = {x = (0,p) , p € T} G T 2 : two points close to 70, A:= (e,p) and 
B:= (1 —s,p), have images that differ, in the e — ► limit, by a vector 
d ( _ 1} (A,B) = (a, a) (mod 1). 

(2) Inverse Sawtooth Maps {S^ 1 } are discontinuous on the subset 

7_i := S a (70) = {x = (p,p) , p G T} G T 2 : two points close to 7_i, namely 
^4:= (p + e,p — s) and i?:= (p — e,p + e), have images that differ, in the 

e — > limit, by a vector d^ii (A, S) = (0, a) (mod 1). 

j a 

(3) The maps T a and T" 1 are continuous: 

a G Z =► (A, B) = dj^ (A, B) = (0, 0) (mod 1). 

(4) The eigenvalues of the matrix S a = ( 1+ a a \ ) are (a + 2 ± ^(a + 2) 2 - 4) /2. 
They are conjugate complex numbers if a G [—4,0], whereas one eigenvalue 
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A > 1 if a [—4,0]. In this case, distances are stretched along the direc- 
tion of the eigenvector |e + ), S a \e + ) = A|e + ), contracted along that of |e_), 
S a \ e —} = A _1 |e_): log A is a (positive) Lyapunov exponent. 
For such a's all periodic points are hyperbolic [18]. 




Figure 2: In the upper row, we depict the effects of the discontinuities of a SM with 
a = i; the picture in the middle shows the discontinuity lines 70 and 7_i, whereas those 
on the right and left show how they evolve backward and forward in time. The different 
parallel bands help the reader to figure out the toral periodicity and the discontinuous 
character of the map, also highlighted by the aperiodic splits of two spots. Further, for 
sake of comparison, the lower row presents the same case of the upper one but for the 
continuous dynamics (a = 1). 



Remarks 4.2 

Because of the presence of the fractional part in (7) and (9), we have to 
distinguish the action of S a and 5" 1 from a mere matrix action. We shall 
adopt the following notations: 

i. With S a the matrix \ ) in Property 4.1.4, the expression S a (x) will 
denote the action represented by (7), whereas S a ■ x will denote the matrix 
action of S a on the vector x. 

ii. When the dynamics arises from the action of the UMG (see Remark 4.1. v.), 
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so, in particular, when {T a } a£Z is the family of toral automorphisms, equa- 
tion (7) assumes the simpler form T a (x) = T a ■ x (mod 1). 

iii. Analogously, expressions like T a ■ x, ■ x, T~ l ■ x and (T^ r ) 1 • x, will 
denote the actions by T a itself, its transposed, its inverse and the inverse 
of the transposed, respectively. 



4.2. Algebraic description of continuous and discretized Sawtooth Maps 

In this Section we make use of the commutative (Von Neumann) algebra (T 2 ) intro- 
duced in Section 2 and consider the algebraic description of Sawtooth Maps by triples 
(Ljf (T 2 ) ,w M ,e a ), where has been defined in (1) and 9 Q : (T 2 ) i-> L~ (T 2 ) is 
the discrete-time dynamics generated as follows: 

e a (f)(x):= f(S a (x)) , a G R • 

The maps @ J a , j £ Z are automorphisms of (T 2 ) and leave the state lo^ invariant. 

Our aim is now to define a suitable discrete evolution ©Af,a on Vj\f, such that the 
discretized triplets (%,Tj\f,9^ a ) converge to the continuous SM. 

We start by introducing two different kinds of maps: the first ones, ut\ j g z, are 
defined on the torus T 2 ([0, A) 2 ), namely [0, N) x [0, N) (mod A"), and given by 



NS° a [^) GT 2 ([0,A) 2 j , (10a) 

r 2 ([o,A") 2N 

T 2 ([0, A) 2 ) B x -C^ (x) := U±\ U±\ ■■■ U^H Ut\ * ))•••)) , 3^ 



T 2 ([0, A) 2 ) Ul{x):= x 



T 2 ( [0, NY)3x ^U^ (x) := N [-) G T (JO, N) j , (10b) 



j times 



= NS*{±) GT 2 ([0,AT) 2 j • (10c) 
The second class consists of maps Va j from T 2 ( [0, A) 2 ) onto its subset (Z/AZ) 2 , whose 



Continuous Limit of Discrete Sawtooth Maps and its Algebraic Framework 15 



actions are as follows 



T 2 ([0, AO 2 ) Bx^ V° (x) := [x\ 

= ± l±U°{[x\)\ G (Z/NZ) 2 , (11a) 
T 2 ([0,N) 2 ) B x^V^ 1 (x):= ± [±U± l ([x\)\ G (Z/NZ) 2 , (lib) 



T 2 



( [0) N f) 3I h ^ ( X ) -.= y^jv+H- • • v*Hv?\ L»J ))•••) ) , 

7 times 

± l± t^ 1 ( ± L± c^ 1 ( • • • ± L± ^ ( ± L± ut 1 ( \x\ ) J )j • • • ) J )j g (z/ivz) 2 • (iic) 



j times 

Remark 4.3 

The maps Ui are extensions of the S J a on the enlarged torus T 2 ^[0, A^) 2 ); 
however, they do not map the lattice Ln into itself, therefore we are forced 
to use the maps V£ to define a consistent discretized dynamics. 

Definition 4.1 

8aT,q will denote the map: 

V X 3X^Q M<a (X):= X Va{l) y a{£) \l){t\eV N ■ (12) 

Qj\f a is a *-automorphism of T>j\r; indeed, the map 

(Z/NZ) 2 B£i — >V a (£)e (Z/NZ) 2 
is a bijection, so that (12) can be rewritten in the more convenient form 

©A/>P0 = E X V a (£),V a (£)\ £ )( £ \ = 



(see Remark 4.4.iii. below) = W t 



E X^lV-'is))^- 1 ^) 

Va 1 {s)e{z/Nz) 2 

( \ 

E Xs,s\ S)(s\ 

Vail equiv. 
classes / 



a.N 



w: N = (13) 



W a , N X W* N , 
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where the operators W a> N, defined by linearly extending the maps 

H M 9 \l) W a>N \t) := | V- 1 {£) ) G H M ■ (14) 
to Htf, are unitary: N \i) := | (-£)). 

For the same reason the state r/y is 0^ ^-invariant and is invertible too. 
Note that Qj^ a := 0_/y a o • ■ ■ o Q is implemented by (•£) given in (11c). 

V '■ v ^ 

j times 

Remarks 4.4 

i. The double ± sign in front and within every floor function in equations (11) 
is needed in order to have Va\ (x)) = V®(x) (the identity when 
x G (Z/iVZ) 2 ); the reason is that, in general, [—x\ ^ — \x\ , for x g" Z 
(see [27]). 

ii. When a G Z, (Z/iVZ) 2 9 i i — ► V Q {£) = T a ■ £ G (Z/NZ) 2 , namely the 
action of the map V a becomes that of a matrix (mod N). Moreover, in 
that case, U a and V a coincide. 

iii. Since £ i — > V a (£) is a bijection, in (13) one can sum over the equivalence 
classes. 



5. Continuous limit of the dynamics 

One of the main issues in the semi-classical analysis is to compare if and how the quantum 
and classical time evolutions mimic each other when a suitable quantization parameter 
goes to zero. 

In this article we are instead considering the possible agreement between the dy- 
namics of continuous classical systems and that of a class of discrete approximants. In 
practice, in our case, we will study the difference 

ei-Jooxoe^o^ (15) 

which represents how much the discrete dynamics at timestep j differs from the contin- 
uous one at the same timestep. 
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For quantum systems, whose classical limit is chaotic, the situation is strikingly 
different from those with regular classical limit. In the former case, classical and quantum 
mechanics agree, that is a difference as in (15) is negligible, only over times j which scale 
logarithmically (and not as a power law) in the quantization parameter. 

As we shall see, such a type of scaling is not exclusively related with non-commutativity; 
in fact, the quantization-like procedure developed so far, exhibits a similar behavior when 
N — > oo and we recover (L^° (T 2 ) , u^Oq.) as a continuous limit of (T>j\f, r/v, @j\f, a ) ■ 

5.1. Continuous limit for Sawtooth Maps 

Later on we shall show that the difference in (15) goes to zero in a suitable topology; 
for the moment we just note that the major difficulties in the proof are due to the 
discontinuous character of the fractional part that appears in (7). 

It is therefore important to briefly discuss the discontinuities of the maps S a [16-18]. 

As already noted in Property 4.1.1, S a is discontinuous on the circle 70; therefore S™ will 
be discontinuous on the preimages 

7m:=S- m ( 7o ) for 0^m<n, (16a) 

whereas the discontinuities of S~ n lie on the sets 

7 _ m :=S a m ( 7o ) for 0<m^n- (16b) 

Apart from 7_i, whose projection on the [0, l) 2 square is its diagonal (see Fig. 5), each 
set of the type 7 m (for 7 _ m the argument is similar) is the (disjoint) union of segments 
parallel to each other whose endpoints lie either on the same segment belonging to 7p , 
p < m, or on two different segments belonging to j p and 7 p /, with p' ^ p < m [17]. 
It proves convenient to introduce the discontinuity set of 5^, 

n-l 

T 2 DT n := (J 7p , (17) 
p=0 

and its complementary set, G n := T 2 \ T n . 

We now enlarge the previous definition from continuous Sawtooth Maps, to discretized 
ones. 
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Definitions 5.1 



We shall call "segment", and denote it by (A,B), the shortest curve joining 
A, B S T 2 , by I {j v ) the length of the curve 7 P and by 



7r 



(e) := \x E T 2 | d T 2 (as, lp ) < e} (18) 



the strip around 7 P of width e, where the distance d T 2 (•, •) on the torus has 
been introduced in Definition 3.1. 



Further, we shall denote by 



n-l 



T n (e):=(j %(e) (19) 

p=0 

the union of the strips up to p = n — 1 and by (e) the subset of points 



G^(e):=^x€T 2 ^0r n (e)j , (20) 

where the lattice points have been introduced in Definition 3.2. 

As already observed, in order to prove that the discretized SM tend to continuous SM 
when N — > oo, the main problem is to control the discontinuities. It proves convenient 
to subdivide the lattice points in a good and a bad set and show that, on the former, 
Va —Ua, at least on a certain time-scale (see Remark 4.3). This will not turn out to be 
true for the bad set, however we shall show that the latter tends with N to a set of zero 
Lebesgue measure and thus becomes ineffective. 

Following this strategy, we shall concretely show that the difference (15) goes to 
zero with N — > oo in the strong topology over the Hilbert space L 2 (T 2 ). More precisely, 
we have the following theorem 

Theorem 1 

Let iDtf, taT) @AT,a) be a sequence of discretized SM as defined in Section 4-' 
for all 7 > 3, 

V/a»(T 2 ) , s-lim (ei-J ooM oQ^ a oJ MiOO )(f) = , (21) 
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where the limit is in the strong topology over the Hilbert space L 2 (T 2 ) and 

rj > y/2 is the largest eigenvalue of the matrix \S a \ : = \J SaS a , with S a 
defined in Property 4.1.4. 

The previous Theorem indicates that the time limit and the continuous limit do 
not commute. In particular, the difference between the discretized dynamics and the 
continuous one can be made small by increasing N, while it becomes large beyond the 
time scale j ~ ^ l " f ^ . This phenomenon is the same as in quantum chaos and points 
to discretization of phase space (in the traditional semi-classical treatment of quantum 
systems), rather than to non-commutativity, as the source of the so-called logarithmic 
breaking time. The constant 7 is a form factor, which reflects the fine structure of the 
dynamics: for instance, in the case of quantum cat maps [25], 7 = 2. 

Remark 5.1 

The parameter 7 > 3 in Theorem 1 may seem overestimated if compared 
with the case of the quantum Cat Map, where 7 = 2. As we shall see (in 
particular in the next Proposition 5.2), the upper bound for 7 is dictated by 
the discontinuities of the Sawtooth Maps, and not by commutativity. The 
corresponding exponent assumes the lower value 7 > 1 in the case of dis- 
cretized Cat Maps, that include Sawtooth Maps with integer a. This result 
will be presented in a forthcoming paper [28], in which we study the breaking 
time tb (N), here ^7^^, relative to the chaotic or non-chaotic properties of 
the dynamics. In particular, in the hyperbolic regime, the parameter log 77 of 
Theorem 1 is replaced by the Lyapunov exponent log A whereas, in the elliptic 

regime, the two limits j, N — ► 00 do commute and in the parabolic one, the 

1 

breaking time is given by tb (N) = N~t . 

The proof of Theorem 1 consists of several steps, among which the most important 
is a property, satisfied by our choice of Lattice States, which we shall call dynamical 
localization. 

We give a full proof that our choice of Lattice States satisfies such property, since it 
represents a natural request that should be fulfilled by any consistent discretization/de- 
discretization (quantization/de-quantization) scheme. 

Remarks 5.2 

(1) In analogy to the quantum case, Dynamical localization is what one ex- 
pects from a good choice of states suited the study of the continuous limit: 
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in fact, it essentially amounts to asking that LS remain decently localized 
around the continuous trajectories while evolving with the corresponding 
discrete evolution. As we shall see this is the case only on logarithmic 
time-scales. Informally, when N — > oo, the quantities 

Kj(x,y) :={C N {x),wl N C N {y)) 

should behave as if M\Kj(x, y)\ 2 ~ S(S J a x — y): this would make the 
discretization analogous to the notion of regular quantization described in 
Section V of [29]. Actually, with our choice of LS, the quantity Kj(x,y) is 
a Kronecker delta. 

(2) In quantum chaos, instead of seeking for the dynamical localization, one 
can study the dynamical spreading of Coherent States. Consider for instance 
the classical function / over the phase space, its corresponding quantum ob- 
servable Opjj (/) and a Coherent State \Cn(x)) centred at the point x. The 
time needed for the quantum mechanical expectation {C^{x), Op/j (/) C%{x)) 
to converge to the average of / over a suitable invariant measure can be 
explicitly analyzed. Recent work [7, 9] shows that also this time scales 
logarithmically in h, at least for the automorphisms on the 2-torus. 

(3) The constraint j < C log M is typical of hyperbolic behavior with Lyapunov 
exponent log A and comes heuristically as follows: the expansion of an initial 
small distance 6 can be exponential until the distance becomes the largest 
possible, namely <5A Tb ~ 1. After discretization, the minimal distance gives 
S = ^, therefore one estimates Tb ^ ^y, which is called breaking time 
and sets the time-scale over which continuous and discretized dynamics 
mimic each other. 

(4) In quantum chaos, the semi-classical analysis leads to an estimate of Tb 
exactly as above; further, the logarithmic dependence on h of Tb is a sig- 
nature of the hyperbolic character of the classical limit. Conversely, if the 
classical limit is regular, then the time scale when quantum and classical 
behaviors are more or less indistinguishable goes as h~ b , b > 0. Another 
interpretation of the breaking time is given in [8], where it is related to the 
shortest time needed for the system to transfer all scales 1 ^ I ^ H down 
to the "quantum scale" h. Indeed, this is the scale at which the differences 
among quantum and classical mechanics come up. Regarding the SM, the 
hyperbolic case corresponds to S a with eigenvalue A > 1, whereas the reg- 
ular cases are the elliptic one (two complex eigenvalues) and the parabolic 
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one (only one eigenvalue = 1). 



(5) The dynamical localization property has fruitfully been used in several 
quantum contexts [25]; however, to our knowledge, this is the first instance, 
though not properly quantal, where dynamical localization is fully exposed. 



Before proceeding with the proof of Theorem 1, it is important to notice that in its 
statement the Lyapunov exponent log A does not appear but log?], instead; of course A 
and 77 are related for A is eigenvalue of S a , and rj of \j SaS a (see Remark 5.1). 

As will become clear during the proof, the use of rj and not of A is required by 
the discontinuous character of SM. In fact, the discontinuities do not allow us to control 
the difference between the n-th iterates of the discretized and the continuous dynamics, 
but instead force us to estimate that difference at each single time-step up to n and 
to put all the estimates together. In the single time-step estimate, independently of 
whether the map is continuous or not, one must use rj, which coincides with A only 
when the dynamical matrix S a is symmetric. Indeed, Figure 3 shows that the eigenvalue 
rj correctly describes how volumes behaves under a single application of the dynamics, 
whereas A underestimates it. On the contrary, it is A™ which asymptotically controls 
the stretching, whereas r? n largely overestimates it. In the regular elliptic case, where 
A = and J] ^ \/2, the use of rj gives the impression of hyperbolic stretching, whereas 
the elliptic motion is confined: from the lower strip in Figure 3 it is apparent that such 
hyperbolicity is spurious. 
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Figure 3: In Plots A, B and C we compare the estimates of the (maximum) stretching 
given by the action of the SM S 1 , 1Q and its temporal iterates Sy 1Q (n ^ 5) given by A, 
respectively rj, on a small ball B® of radius v, centered in (|, |) G T 2 . The five evolved 
images of the ball, namely {B™ , 1 < n ^ 5}, are plotted together with B®, using different 
colors. In A we surround every evolved ball B™ with the smallest circle containing it. We 
compare that plot with B and C, in which the surrounding circles have radii proportional 
to X n v, respectively rj n v; in both cases the correct radii of A are overestimated although, 
on the long run, circles in B provide a good approximation. 

The fake hyperbolicity given by r\ is clearly shown in D and E, where a parabolic SM So 
and an elliptic one SLi^o are presented: in the first case the maximum spreading grows 
linearly, whereas in the second one it remains confined, and the estimate given by the 
surrounding circles of radii growing as powers of rj is inappropriate. 
Note that in all examples C-E, the black circles of radii rj v rightly surround B^. 



Theorem 2 (Dynamical localization with {\Cj^(x))} states) 

For a G R, G K+ \ (0, 2 ] and d > 0, there exists N = N (a, /?, d ) G 1N+ 
with the following property: if N > No and n < ^ , then 

d T 2 (52 (x) ,y)^d ^ (C M (x) | W2 )N C K {y)) = , 
for all y£T 2 and x G ( ^ ) , where W™ N is the unitary operator defined 
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m 



(14), N = 2 V2 (V2 + l) rj 2n and G% (e) has been introduced in Defi 
Hons 5.1. 



ni- 



In order to prove Theorem 2, we need the following result, whose proof can be found in 
Appendix B. 



Proposition 5.1 

With the notation of Definitions 3.1 and 5.1, and with [E]° denoting the 
complement of E C T 2 , [E]° := T 2 \ E, the following inclusions hold: 



C G% (e) C 



\/2iV 

Further, for a € K, and n € fJ + , if 



r„ e 



V2N 



(22) 



iV > N = 2 ( V2 + 1) r] 2n and x £ then 



(23) 



Proof of Theorem 2 : 



Using the definition of {\Cj\f(x))} in (6), we easily compute 

(C U (x) | W£ )JV C M {y)) = (x N | I/"™ (y N ) ) = 5 
Using the triangular inequality, we get: 



(TV) 



TV 



(24) 



iV 

, , ( Vg{x N ) Vn\,, (VN 



or equivalently, using the Definitions (10), 
'V?(£jv) Vn 



d^i^f^,^-) ><h*(SZ(x) ,y) 



, US(Nx) V£(x N ) \ (y N N 



24 F. Benatti and V. Cappellini 



Now, since d V 2 (S™ (x) , y) ^ cfo by hypothesis, using (40) in Appendix B and observing 
that x £ ySf) P erm its us to use (23) in Proposition 5.1, namely that 

N>N => d ^{-^^^^^)^w{-^r) ' (25) 

we can derive 

The r.h.s. of the previous inequality can always be made strictly larger than 

( vS{xn) y N \ l , . 

^{-^—'ir) > n ' (26) 

by choosing an N larger than 

N M (n) = max ij- 1 + ^2 ( ?? " ? + ^~ ~ ) + , A = 2 (v^ + l) t? 2 " } , (27) 

so that the condition on the l.h.s. of (25) is also satisfied. From (24) and (26), we have 
N>N u (n) =}► <CV(a;)|W^ jAr CV( 2 /)>=0- (28) 

Indeed, if the toral distance between two points (z, w) exceeds then the corresponding 
grid points (zn,wn) are different and then the periodic Kronecker delta in (24) vanishes. 

Since the (non-decreasing) function Am in (27) is eventually bounded by rf n ((3 being 
strictly greater than two), we define n as the time when Am (n) = rf™ —. Ao, and choose 
A > Ao, x £ G% (Mr) ■ Thus > if < n < % then A' • A n A M (n) > A M (n), whereas 



n \2N ^ 

if n < n < ; then A > rf n > A M (n) and (28) holds for all < n < Jgigf • ■ 

In order to proceed with the proof of Theorem 1, we need another auxiliary result which 
is proved in Appendix C. 

Proposition 5.2 

With the notation of Definition 5.1, the following relations hold for all p G M, 
n E JNT + and e S E+: 

I ( lp ) < 77 p , (29a) 
A* (%(£)) ^2e^ + 7re 2 , (29b) 
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A*(r n (e)) < 2{y/2 + l) erP + irne 



(29c) 



Moreover, if N G JNT+ and N = 2 ^2 (\/2 + l) ?] 2n (cfr. equation (23) in 
Proposition 5.1): 



N > N 



1 1 



G 



N 



N 
2N 



38 rj 
N 



3n 



(29d) 



We are finally in position to conclude with 



Proof of Theorem 1: 



We subdivide the proof in two steps: in the first we concentrate on continuous /, that is 
/ G C° (T 2 ) (c -L 2 (T 2 )); in the second one we extend the result to essentially bounded 
function by applying the following Corollary of Lusin's Theorem [21,30,31]: 

Given f G (X), with X compact, there exists a sequence {f n } of contin- 
uous functions on X such that |/ n | < ||/||oo and converging to f \i - almost 
everywhere. 



(1) Let / G C°(T 2 ) and Op jVV (/):= [&> a - ° ©V,a o J Koo j (/): notice that 

Opj N (/) is a multiplication operator on L 2 (T 2 ), but also an (T 2 ) (and thus also 
an L 2 (T 2 )) function. According to (21), we must show that 

V<?GL 2 (T 2 ) , lim || HN (f) g\\ 2 = ■ 



3< 



1 log JV 

7 log); 



Using Schwartz's inequality first with g in the class of simple functions and then using 
their density in L 2 (T 2 ), we have just to show that 



lim 

j,N—>oo 

■ < 1 log N 
<* 7 log 7] 



Op,-v(/) 



Explicitly, using (1), we write: 

2 a i/;' ()i',.v (/)) 

+ w M 



Op,-v(/) 



2 =UJ V 



(©i/)*(©i/) 



+ 



Joo,AT ° Q ° JV.ooJ (/)* ( Joo,M ° Q o JV^J (/) 

2 Re [(e> /)* (joo,^ o o J N)00 ) (/)] } , 



+ 
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which, via Proposition 3.2.1, becomes 



: ^ H (/) ©£ (/)] " 2 ^ {w [^r,oo (© J Q /)* (©^ J*r,°o) (/)] } + 

+ W {JM,oo° Joo,M oQ Af, a ° ^,oo) (/)* (©A/> ° ^Zv,oo) (/) , 



that, using Proposition 3.2.3, can be recast as 



o e> ) (//) + w o ^ i00 (/)* e^ )Q o ^ i00 (/) 



+ 



2 Re i Tj\f 



{ 



(/)* (©^a ^V,») (/)] } 



^(l/l 2 ) + (wo©irJ [^V-,oo(/r^,oo(/)]-2Re(/^(/)) , 



with 



i,7V (/) := W [(^ |00 o 6i) (/)* (®if,a J*,°o) (/) 



Jt 2 Jt 2 



Now, Proposition 3.2.2 yields 



T M ° ©AT a ) [JV,oo (/)* JV,oo (/)] = W [JV,oo (/)* JV,oo (/)] — ► <*M I / I 

so that the strategy is to prove that also Ij t N (/) goes to uj^ (\ f\ 2 ^j = J fj,(dx)\f( 
when j, N -> oo with j < li&K. 

J1 J -y jog 77 

Resorting to m Definition 5.1, and to its complementary set 



r N ( JL\ 



1,2 \ G n ( m ) ' we can write 



ij, N {f)- / Mdy)|/(y)h 

T 2 



M(dx) / ^dy)f(y){f(Six)-f(y))Af\{C N (x),Wi N C^y))\' 
t 2 Jt 2 



< 



+ 



\ 2JV , 



//(da) / ^dy)f{y){f{Slx)-f{y))M\{C N {x),Wi N C N {y))^ 



T 2 



u n I IN 



f,(dx) / M(dy)/(y)(/(^)-/(y))^KC^(aj),W-^cv(y))p 

Vt 2 



(30) 
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For the first integral in the r.h.s. of the previous expression we have: 



G N 

n I 2JV 



fi(dx) / yi{dy)f{y){f{Six)-f{y))N\{C M {x),Wi N C N {y))Y 



T- 



< 2 



< 2 



C'N N ) 
l 2N ) 



H(dx) / Kdy)M\({W* tN yC M (x),C M (y))\ 2 
Jt 2 



\2 



A' 



2iV 



76^' 



where we have used completeness and normalization Properties 3.1 and equation (29d) 
from Proposition 5.2; this term becomes negligible for large N > N iff j < }■ , with 
7 > 3. 

Now it remains to prove that the second term in (30) is also negligible for large TV: 
selecting a ball B(SaX,do), one derives 



G N -2- 

n \ 2N 



Mdx) / ix{&y)f{y){f{Slx)-f{y))N\(C N {x),W 3 aN C N {y))t 
Jt 2 



< 



+ 



QN JL 
« \ 2N 



G N -2- 
n I, 2N 



»(dx) / . K&V) f(y)(f(S j a x) - f{y))N\{C M {x),Wi N C M {y))\ 2 

JB{S 3 a xA)) 



fx(dx) I v(dy)f(y){f(Six) - f{y))N\(C N {x),W ] aN C N {y))f 

lT 2 \B(S J a x,d ) 



Applying the mean value theorem in the first double integral, we get that 3c £ B(SaX, do) 
such that 



rN JL 

I 2N i 



< 



r N ( JL 



/^(dy) AA|((W"* iJV ) J CV(a;), Cj K f(y))\ 2 



H(dx) /(c) (/(5jaj) - /(c)) 



+ 2 



I 2JV 



M(da;) / n(dy)/V\(C^/(x), W J a N Cj^(y)}\ 2 

lT 2 \B(Si X ,do) 
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Finally, using completeness and normalization (Properties 3.1), we arrive at the upper 
bound 



< 



sup 

ceB(z,d ) 



|(/(*) -/(c))| +2 



2 



sup 



\{C N {xlWi N C N {y))\ 



ygB{S 3 aX ,do) 



By uniform continuity, the first term can be made arbitrarily small, provided we choose 
do small enough. For the second integral, we use Theorem 2, which provides us with 
No = No (do) depending on the same do , such that the second term vanishes for all 
N>N andfarallj<iM 

(2) In order to extend the result of point (1) to / S (T 2 ), we use the Corollary of 
Lusin's Theorem, choose a sequence {f n } n as in its statement and estimate 

. lim [J Vj>N (/) || < . lim || Op jN (f - f n ) || 2 + . lim || Op j N (/„) || • 

1 log N ■ , 1 log N \ log JV 



3< 



7 log V 



J< 



7 log V 



J< 



7 log V 



Using point (1), the second term in the r.h.s. of the previous equation can be 
bounded by arbitrarily small s, indeed f n S C° (T 2 ). 

For the first term we proceed as follows: using Definition 4.1 together with equa- 
tions (37) and (38) of Appendix A, we find 

(JooM ° ° JM,oo) ( 9 )(x) = Yl T " (9) (^4r) X Q N (j,)( x ) > ( 31 ) 



where g is any measurable function on T 2 . Then, because of how the running average 
operator (RAO) T^v is defined, for all g G (T 2 ) it follows that 



[Joo,M ° ©at.q ° (9) 1 ^ {joojj o ® j Ma o Jm )CO ) (I g 



\9\\i 



where || - 1| ^ denotes the (T 2 )-norm, and that 



sup 



JooM 0< S>Af a ° JN,oo (9) 



r^v (g) 



N 



| < lir^^Ho^ \\g\ 



Indeed, the first equality in the last formula comes from the definition of essential 
norm [21] (which in this case amounts to the greater absolute value assumed by the 
simple function Joo,M ° ®j\f a ° «^V,°o)) whereas the first inequality is a consequence of 
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the continuity of Tjq and the last one from Proposition A.l. Putting last two inequality 
together, we obtain 

[Jaajf ° ©AT,a ° ^,00 J (fiO < Halloo llsll 1 , 
whence, setting g = f - f n , 

!! Op,- N (/ - /„) || 2 = I 6i (/ - /„) - JooM ©foe ^,00 (/ - f n ) 1 2 (32) 

< II / - fn || 2 + II / - fn Hoc || / - /n 111 , Vj, N ■ 



Now convergence follows from Lusin's Corollary. 
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Figure 4: These two plots show how the difference between Ja/\oo ° ©a and Oa/> Jm,od 
becomes smaller with N. For the continuous SM, Gi, the actions JV,oo ° ®i and Bjv\i ° 
JN, oo on / G L~ (T 2 ) (left part of both plots) are plotted for two different N: N = 16 
(top) and N = 48 (bottom). The resulting matrices are mapped back, together with the 
function 0i (/), on the unfolded torus, by means of the de-discretization operator Joaj^I- 



Continuous Limit of Discrete Sawtooth Maps and its Algebraic Framework 31 




7_i = 53/2(70): 7-2 = 5" 3 2 /2 (70): 

70: Discontinuity Discontinuity line Discontinuity line 




Figure 5: Here, the same picture as in Figure 4, is represented, with a finer discretization 
given by N = 120 and a different function g 6 (T 2 ), for a discontinuous SM, 63/2, 
acting two times. Choosing a function g with sharp variation across 70 (blue lines), 
the preimage of 7_i, the discontinuity of ©3/2 makes it evident how the differences 
between J7i44oo,oo © Q and ©14400,0 «^i4400,oo are the greater the closer they are to 
the discontinuity line 7_i (red lines). Of course, the longer the temporal evolution, the 
worst the correspondence, in the sense that several new discontinuity lines come to play 
a role. In the case at hands, the map acts twice, and 7-2 is felt by ©144000 i7i4400,oo, 
as expected. 
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6. Conclusions 



In this article we have considered discrete approximants of Sawtooth Maps on the torus 
and we have studied them in an algebraic framework modeled on the so-called Anti-Wick 
quantization; In fact, finite-dimensional discretization and quantization can be seen as 
similar procedures in that they map an abelian Von Neumann algebra (of essentially 
bounded functions on phase-space) into finite-dimensional matrix subalgebras, the only 
difference being whether the latter are diagonal (commutative) or not. 

In the semi-classical analysis of classically chaotic quantum systems, the correspon- 
dence classical/quantum is usually observed only on time-scales that are logarithmic in 
the quantization parameter h. The motivation of our study was to show that the same 
phenomenon arises when a hyperbolic classical system is discretized, namely forced to 
move on a lattice, and afterwards the lattice spacing is sent to zero. 

Previous results [14] based on the numerical investigation of the entropy production, 
indicate that it should indeed be so; however, these results were not supported by a solid 
framework where to analyze the continuous limit of the family of discrete approximants. 
This is the content of this article. 

The major difficulty was represented by the need of controlling the discontinuous 
character of Sawtooth Maps, which was made possible by an appropriate choice of Lattice 
States. In fact, similarly to the entropic approach which, despite the dynamics being 
classical, was based on a quantum dynamical entropy, the discretization/de-discretization 
procedure we set up is based on quantum tools. 

The choice of Lattice States was naturally pointed to by the lattice structure of 
the discrete phase-space and turned out to posses the right localization properties for 
mastering the discontinuities. The result is the appearance of a logarithmic time-scale 
when the discrete hyperbolic SM tend to their continuous limit; namely, the continuous 
and discrete dynamics agree up to a breaking time which is proportional to the logarithm 
of the lattice spacing. 

The proportionality constant does not involve the Lyapunov exponent, that is the 
eigenvalue A > 1 of the dynamical matrix S a , rather the largest eigenvalue, rj, of \J SaS a . 
In the case of elliptic SM, |A| = 1, rj > y/2; however the resulting breaking time is a 
spurious effect, while when A > 1, the presence of rj in the breaking time seems to be an 
unavoidable consequence of the discontinuous dynamics. 
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A. Anti Wick discretization of (T 2 ) 

In this appendix we will apply Definitions 3.3 and discretize (T 2 ) by means of the 
LS set {\Cj\f(x)) | x G T 2 } G Tij^f introduced in Section 3.1. 

In this framework, the discretizing/de-discretizing operators of Definitions 3.3 read: 

L™(T 2 )3f^N 2 [ ^dx)f(x)\x N ){x N \ =:J* l00 (f) eV N , (33) 
V m bX^(x n \X\x n )=:J ooM (X)(x)£S{T 2 )cL™(T 2 ) , (34) 

where S (T 2 ) denotes the set of simple functions [21] on the torus. The matrix elements 
of JV,oo(/) are as follows: 

M i,m : = <* I •&/>,(/) I m) 

= iV 2 / M (dx) /(x) (£\x N ){x N \m) 

= N-jf 1 ^ /'d, 2 /(x) C..„, 

- *• c e, c L ^j ■ 

This implies 

<™ = ^C L-\ ^ L] d ^/W , (35) 

so that varying / G U£ (T 2 ) yields Ran(JXf,oo) = T^M ■ In order to recast (35) into a 
nicer expression, we introduce 

Definition A.l (Running Average Operator (RAO)) 

Let Qm (x) denote the square of side 1/N, oriented parallel to the axis of the 
torus and centered around x; then, the Running Average Operator 
T N : (X) i — ► C° (T 2 ), is defined by 

(T 2 ) 3 f(x) T N (/) (x) = : N 2 [ M (dy) f(y) G C° (T 2 ) • 

JQ N (x) 
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Proposition A.l 



Given / e L£° (T 2 ) , the function :=T N (/) is uniformly continuous on T 2 ; 
moreover, the Running Average Operator has norm 



|IV|| B := sup 



\Tn (/)llr 



1 • 



(36) 



Proof: 



Let xq £ T 2 , x S Qn{xo) and Xe denote the characteristic function of E C T 2 . By 
Definition (A.l): 



(-o) - f™ (x 



N 



H(dy) f(y) {X QN(xQ) (y) - X QN{x) (y)) 
o / Mdy) \X QN{xo) (y) - X QN(x) (y)\ 



N 



>N (X 



, ^(Qn (x ) U Q n (x)J - n\Q N WnQ 
According to our hypothesis, x S Qat (^o), thus geometrical considerations lead to: 
h(Qn (xq) UQ n (x)^J < + | xi - x i | ) (— + | x 2 - £ 2 



A*(Q;v (a5o) n Qat 
AmQat (so) U Qn (x) 



— -IXi-Xoil^— -|X 2 -X 02 | 



/■*( <2tv (aso) n Qat (as) ) ^ — ( | xi - x i | + | x 2 - x 2 



2^ 



\x - x\\ , 



so that 

rW) 



< 2y/2 N 



\xq — x\ \ , which proves the continuity of 



f N , while uniform continuity comes from T being compact. 

Concerning the norm in (36), the upper bound ||rjv|| B ^ 1 is clear and the maxi- 
mum is reached by choosing / constant. ■ 



By means of the RAO, the discretization operator in (33) can be conveniently written as 

«Zv,oo(/) = E fP(jf)\')(*\- ( 37 ) 
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Analogously, the de-discretization operator in (34) can be recast as 



J 00 , M (X)(x) = £ X t4 4%= £ X £>£ X QN(i) (x) , (38) 



£e(z/7VZ) 2 £e(z/JVZ) 2 

thus proving that Ran (J^oojv) = $ \ ^ 2 ) ■ 

Moreover, combining equations (37) and (38), we explicitly get the simple function 
arising from / G (T 2 ), via AW discretization/de-discretization: 

(Joo,M°JM,oo)(f)(x)= Yl TN ^(jf) X Q»(ir)W- (39) 

ie(z/NZ) 2 ^ ' 

The action of the operator Joo,N Jm,oo can be seen in Figures 4 and 5. 



B. Proof of Proposition 5.1 

We start by proving the inclusions (22). 

For every real number t, we have ^ (Nt + \ ) = Nt + \ — \ Nt + < 1, so that 

[Nt + \\ 



N 



From (5) in Definition 3.2, we derive 



2N 



d V 2 x 



x N 



1 



' N J y/2N 
Then, let us consider the triangular inequality 



d T 2 (x , y) ^ d T 2 ( x , ^ 



<*T 2 ( -^r , y 



and let us take the infimum over the set y G T n defined in (17) 



(40) 



(41) 



V ( ^ , r n j > d T a (x , r n ) - d T 2 ( x , ^ 



^ d T a (x , r n ) 



i 



V2N 
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where we used (40). Therefore, considering the complement [r n (e)]° of the union of 
strip of width e, T n (e) defined in (19), we get that 



* e[ r„(C =* *s 6 [r„(.--^) 



Further, from (20), it follows that, if the lattice point ^ does not belong to T n ye — 

then the corresponding point x G T 2 must belong to [e — ■ 

Changing e— i — > e we obtain the first inclusion relation in equation (22); the second 

one follows by interchanging the role played by ^ and x in (41). 

In order to prove (23) , we start by considering the matrices S a = ( 1 + a \ ) and 



its inverse S^. 1 = (j^i+q)- Let r\ be the largest (positive) eigenvalue of y SaS a ; its 
characteristic polynomial for rj is if — (2 a 2 + 2 a + 3) rj 2 + 1 = 0, whence rj attains its 
minimum 77 m ; n = y/2 at a = —\. Then, we set N:= 2y/2 (\/2 + l) rj 2n , n S M, choose 
N > N and proceed by induction. 

p = : from definitions (10) and (11), it follows 



"^v— '^v— ) =d ^ { X >1T) < V^ < JT 



where the first inequality follows from (40), thus relation (23) holds for p = 0. 
p = q — 1 , l^q^n: since 

* iV ' — }V — J ^ ^ ( n ' N 1 



+ d V 2 



N ' N 



using (10) in the first term and noting that, from definitions (10) and (11), the second 

/2 

term is less or equal to -j^, we get 

f U q a {Nx) V«(x N ) \ < / f ut 1 (Nx) \ s ( VdT 1 {x N )X\ , V2 



N ' N J * \ \ N \ N M N 
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By the induction hypothesis we have: 

d ^ 2 { — n — ' — n — j^^V^rr) (42) 

<r ^ 1 g 
7 = n H 

f]> V2 , l^q^n ) < - rf~ 2n < -rf 1 ■ (43) 

Now we set e = taking into account that rj ^ s/2 and use the right inclusion in (22) 
to deduce that 



2NJ *~ " V 27V y^iV y ' 

At this point, we make use of the following result, which shall be proved in Lemma B.1.3: 
it states that if a point does not belong to T n (e), the union of the the strips of width 
e ^ \ up to time n, then its orbit under S a up to time n — 1 is farther away than 
£7]~ q , ^ q < n from the discontinuity line 70. Explicitly 



x£T n (e) =► d T 2 (S q a (x) , 7o ) > e rf q , V < q < n , 



whence 



TV V / 

>£(£r). 

r? n — 1 1 

where the second inequality comes from 77 ^ v2\ the relation ^ —= if 1 and 

J? - 1 V2 - 1 
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the following estimates: 



v^iV 



V2 



N 
N 



V 



2 + 1 

2n-l 



V2-1 



- 1 1 

- + 2 



„2n 



77 + 77 — 1 



V2(V_ 



7] — 1 
2n-l 



1 

2 

-.9-1 



7] — 1 



Therefore, comparing (44) with (42) 

^ iV ' iV J <d ^ jj 'To J , Vg^n- 

As a consequence, the segment (^ Ua ^ Nx \ Va j^ Xn ^ cannot cross the line 70. This 
condition, together with (43), allows us to use another result proved in Lemma B.l.lb, 
which states that if a segment (A, B) on the torus does not cross the discontinuity line 
70 then d T 2 (S a (A) , S a (B)) ^ 77 d T 2 (A, B). We can finally conclude with: 

N ' N ) ^ V N \ri-l J N N \ 7]-l J ' 

The following Lemma, which has been used in the proof of the previous Proposition, 
deals with the geometrical properties of the Sawtooth dynamics. 




Lemma B.l 



With 77 the largest (positive) eigenvalue of y S a S a and A, B G T 2 such that 
d T 2 (A,B) < \ r/ -1 , it follows: 

(la) If the segment (A,B) does not cross 7_i, then 

d T2 {S- 1 (A) , 5" 1 (B)) < 77 d T 2 (A, B) ■ (45a) 

(lb) If (A, B) does not cross 70, then 



d T 2 (S a (A),S a (B)) < v d T 2 (A,B) 



(45b) 
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(2) For any given a6K,p£ M + and ^ e ^ \ rj 1 , 



x 



G lp-i (e) => S a 1 (a;) G (t p (t) e) U 7o (r? e)) 



(3) For any given a £ E, n £ M + and ^ e ^ |, with as in (10), 

'Ul{Nx) 



x(£T n (e) =>• d T 2 



N 



Proof: 



In the course of the proof, we shall use that 



l^ 1 -«|| R2 > ^IMb? , 



(46a) 
(46b) 



which directly follows from the definition of 77, where v is any 2-dimensional real vector. 

In order to prove (45), it is convenient to unfold T 2 and the discontinuity of S a 
on the plane K 2 . This is most easily done as follows. Points A G T 2 = ]R 2 /Z 2 are 
represented by equivalence classes 



[a}:= {a + n , n G Z 2 } , a G [0, 1) 



(47) 



Given A,BgT 2 , let A b G [a] be such that 



Notice that 



d T2 ([«],[&]) 



d T 2 ([a] , [b]) = \\a - 6| 



A°-b 



R 2 



R 2 



iff ||a — b|| R2 ^ 



(48) 



(la) (-A, S) not crossing 7_i means that the segment (A b ,fy does not intersect 7_i. 
Periodically covering the plane-K, 2 by squares [0, l) 2 , the 7_i-lines form a set of (parallel) 
straight lines x\ — X2 = n G Z; it follows that b) does not cross 7_i iff 



4-4 = L6i-6 2 J , 



(49) 



where the integral part on the r.h.s. takes values 0,-1, depending on which side of the 
diagonal 7_i the point b lies within. 
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As are not sensitive to the integer part of their arguments, their actions are the same 
on all elements of the equivalence classes (47) , that is 



d T 2 (S- 1 (A) , S" 1 (B)) = d T2 {S- 1 ([a]) , S- 1 ([&])) = d T 2 fc 1 U b ) , S^ 1 (b) 



By expanding (x) = x — [x\, using the definition of S& 1 (•) and putting together all 
integral contributions, condition (49) yields 



d T 2 {S~ l (A) , S" 1 (B)) = min || S~ l (A) - S" 1 (B) + m || r2 

j S- 1 ■ [A b -bj+m' 
d T 2(S- 1 -(A b -b),0) • 



min 



R 2 



Applying (46a) , since we assumed d T 2 (A, B) < ^ rj 1 , we estimate 



Si 1 - A 



R 2 

77 d T 2 {A, B) < - ■ 



A b -b 



R 2 



In particular, using (48), the previous inequalities imply 



d T2 (sz 1 -{A b -b),0 



SI 1 - A 



R 2 



< r? d T 2 (A, 5) 



(lb) Using the same argument as (la), the union of 70-lines constitute a set of straight 
lines x\ = n £ Z; Therefore the segment (A 6 , &) does not cross 70 iff 



4| = L^iJ • 

As done before, by means of (50), we arrive at 

d T 2 (S a (A) , S a (B)) = d T 2 (S a (A b ) , S a (b)) = d T2 (s a -(A b -b),0 



(50) 



The proof can now be completed exactly as for point (2a) before. 

(2) We denote by d V 2 (a;, 7) = inf d V 2 (x,y) the distance of the point x G T 2 from a 
curve 7 S T 2 . Then, from Definition (18) we have: 



x G 7p-i ( e ) => e > d T 2 (aJ, 7p-l) = d T 2 2/" 



(51) 
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where y* is the nearest point to x belonging to 7 p -i. 
We distinguish two cases: 

(2') The segment (x,y*) does not cross 1 7-1 

From (51) and point (la), since S~ (y*) G 7 P (see (16a)), we get 

d T2 {S- 1 (x) , lp ) d T * {S- 1 (x) , S- 1 (y*)) 
< ?? d T 2 (x, y*) < 77 £ • 

Therefore S^ 1 (a;) G 7 p (77 e)- 

(2") The segment (x,y*) crosses 7_i. 

In this case, there exists z G 7_i such that 

d T 2 (x, y*) = d T 2 (x, z) + d T 2 (z, y*) • (52) 

Then, from (51) and (52), 

e > d T 2 (x, y*) ^ d T 2 (x, z) • 
Since, according to (16), S" 1 (z) G 70, from point (la) we get 
d T 2 (S- 1 (x) , 70) < d T 2 {S- 1 (x) , 5" 1 (z)) < r? e , 
that is 5 1 " 1 (x) G 7 (??e). 
(3) From point (2), it follows that, when ^ e ^ \, for p G M + , 

x (% (e) U 7 (e)) 5 Q (x) 7^ (t^e) • (53) 
We prove by induction that, when < e ^ ^, for m G M + , 

m m— 1 

* ^ U % (e) => 5 Q (x) [J 7 P fa -1 *) • (54) 

p=0 p=0 
For m = 1, (54) follows from (53); if (54) holds for m = r, then take 

r+l r 

x 7 p (e) . This means that x ^ |J 7 p (e) and x G" (7 r+ i (e) U 7 (e)) ■ 

p=0 p=0 



1 we stipulate that, if y* € 7-1 or x 6 7-1, we are still in a non-crossing condition 
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Now, using the induction hypothesis and (53), we get 

r+l ) — 1 

x & U % (e) => (a:) |J % {v" 1 ?) and S a (as) % (rj^e) ■ 

p=0 p=0 

Setting m = n — 1 and iterating q times the implication (54) argument, we get 

n— 1 n— 1— q 

a^(J 7p(£)=><S«(*)2 U %( ? r^) , V0<g<n- 

In particular 5^ (a;) 7 (?? _<7 e) , which leads to the lower bound 

d T 2{S q a (x), l0 ) >r]- q e , V0^g<n, 
whence the result follows in view of Definitions (10) and (19). ■ 



C. Proof of Proposition 5.2 



(a) In (16a), we have defined 7 P = Sa P (70) where S~ l (x) (as well as Sa P (x)) is a 
piecewise continuous mapping onto T 2 with jump-discontinuities across the 7 p lines due 
to the presence of the function (•) in (9). Away from the discontinuities, S a p (x) be- 
haves as the matrix action Sa P • x. We want now to estimate the length l(jp); in 
order to do that, we unfold 7 P on the plane and calculate the length of the segment 
| a; € H 2 x = Sa P ■ (y) , y G [0, 1) j, which, in its turn, is the image of 70 under the 
matrix action given by Sa P • x. Therefore, using (46a), the result follows. 

(b) Let L (e) denote the set of points having distance from a segment of length L smaller 
or equal than e: it has a volume (under the Lebesgue measure fi) given by 



2 



where the last term on the r.h.s. takes into account rounding of the extremes of the strip 
by to semi-circle of radius e. Then (29b) follows from (29a). 
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(c) This follows from Definition (19): 



'n-l 



n-1 



v (r n (e)) = ^(U%( e )HE' i (% 

,p=0 / p=0 



Using (29b), we can write: 



n— 1 n— 1 



/' 



(r„(e)) ^2e^? ? p + ^7re 2 = 2e-^— y+nvre 



p=0 p=0 



Finally, the estimate ^ ( v2 + 1 ) x p , valid for x > y/2 yields 

x — 1 V / 



/i(T n (e)) <2e(V2 + l) ?f + nvr£ 2 



(d) By writing the left inclusion in (22) in terms of complementary sets, with e 
we get: 



N 

2N ' 



r N (£_ 

n \ 2N 



- N 1 

c r„ I — + 



G 



iV 



N 
2N 



2N ^2N 



and so 



AT 1 
+ 



2N ^/2N 



By substituting in (29c) N 2?f* = jff + 7/|]y m the place of e, we get: 



G " 1 2iV 



iV 



n 



7T ■ 



V2 + 1 2iV 



Finally, the r.h.s of (55), can be estimated by the following upper bounds: 

N + V2 



7T ■ 



2N 
n 



< 2 



< 77 " 



v^+l 
iV + v^) (v / 2 + l) < 19 



,2'( 



which hold for V N > A 7 ", 77 ^ \/2 and V n E M + . This ends the proof. 



(55) 
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